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ABSTRACT 


The  fundamental  goal  of  this  thesis  is  to  determine  the  geoacoustic  parameters  of  a 
shallow  water  seabed  using  direct  analysis  methods  on  ambient  noise  and  broadband 
explosive  acoustic  data.  All  data  considered  are  from  the  Mid- Atlantic  Bight  shelf  break 
experiment  that  was  conducted  from  19  July  to  9  August  1996.  Simple,  theoretical 
treatments  of  acoustic  propagation  in  a  shallow-water  waveguide  are  applied  to  specific, 
measurable  quantities  in  the  data  which  can  be  inverted  directly  to  produce  estimates  of 
bottom  compressional  sound  speed,  density,  and  attenuation.  Shear  influences  are 
neglected  throughout.  Specifically,  vertical  coherence  of  the  ambient  noise  is  used  to 
determine  the  soimd  speed  contrast  at  the  water/bottom  interface,  mode  travel  times 
extracted  from  spectrograms  of  explosive  data  are  used  to  estimate  bottom  density  based  on 
the  concept  of  an  ideal  waveguide  effective  depth,  and  mode  attenuation  as  a  function  of 
range  extracted  from  similar  spectrograms  are  employed  to  estimate  attenuation.  These 
direct  inversion  methods  are  less  accurate  than  sophisticated  matched  field  processing 
techniques  or  direct  core  measurements,  but  they  do  provide  a  relatively  simple  means  of 
obtaining  reasonable  estimates  of  ocean  bottom  parameters  from  minimal  information. 
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I.  INTRODUCTION 


A.  MID-ATLANTIC  BIGHT  OBJECTIVES 

Better  knowledge  of  operating  environments  may  allow  any  Navy  to  improve  the 
use  of  its  USW  systems  in  different  tactical  scenarios.  Even  though  undersea  soimd  and  its 
practical  uses  have  been  continuously  explored  since  World  War  H,  testing  of  other 
techniques  are  necessary  in  order  to  obtain  important  sea  environment  information  which 
may  help  in  die  development  of  accurate  Naval  USW  systems.  To  face  this  challenge,  the 
Mid-Atlantic  Bight  (MAB)  experiment  was  conducted.  The  MAB  experiment  was  an 
intensive  field  study  of  the  shelfbreak  south  of  Long  Island  during  the  strongly  stratified 
summer  season  whose  basic  objectives  were: 

1.  To  provide  detailed  time  series  of  the  oceanographic  conditions  along  typical 
major  acoustic  paths  by  using  moored  acoustical  arrays. 

2.  To  provide  a  diversity  of  frequencies  and  source  locations  from  broadband  SUS 
sources. 

3.  To  deteimine  the  temperature  structure  of  the  ocean  environment  by  using 
acoustic  tomography. 

This  study  may  provide  the  Navy  an  alternative  and  simple  means  of  determining 
characteristics  in  unknown  or  poorly  known  operating  environments.  (Pickart  et  al.,  1996) 

B.  STUDY  OBJECTIVES 

The  objectives  for  reaching  the  fundamental  goal  of  this  study  are  as  follows: 

1.  To  determine  the  bottom  sound  speed  by  measuring  the  coherence  of  two 
vertically  displaced  sensors  in  the  shallow  water  channel  and  comparing  this 
function  with  one  obtained  from  a  simple  model  equation  over  a  range  of 
frequencies. 

2.  To  determine  the  bottom  density  by  comparing  mode  travel  times  of  computed 
spectrograms  for  a  phone  from  SUS  explosive  data  and  those  obtained  by 
employing  the  concept  of  an  effective  depth  for  an  ideal  waveguide. 
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3.  To  determine  the  bottom  attenuation  by  finding  mode  attenuation  as  a  function 
of  range  extracted  from  spectrograms  and  relating  it  to  loss  per  bottom  bounce. 
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n.  EXPERIMENTAL/ENVIRONMENTAL  DESCRIPTION 


The  summer  Mid- Atlantic  Bight  experiment  consisted  of  a  cruise  from  19  July  to  9 
August  1996  aboard  the  RAA  Endeavor  from  the  University  of  Rhode  Island.  The  cruise 
included  three  major  components:  acoustic  tomography,  SEASOAR  operations,  and 
hydrography/tracers.  An  AXBT  flight  was  also  conducted  during  this  time  period,  as  well 
as  a  SUS  drop  fly-over.  The  acoustic  portion  of  the  experiment  consisted  of  three 
components: 

1 .  The  moored  acoustic  transmission  array,  or  tomography  array 

2.  The  moored  physical  oceanogr^hy  array,  which  provided  detailed  time  series 
of  the  oceanogrq)hic  conditions  along  our  major  acoustic  path 

3.  The  broadband,  explosive  (MK61  SUS  shot)  source  component. 

Multiple  lines  of  SUS  drops  were  deployed  in  the  area  by  P-3  aircraft  and  monitored  by  two 
vertical  line  arrays  (VLAs)  moored  on  the  shelf  in  ~85m  water.  (Pickart  et  al.,  1996) 

To  achieve  the  objectives  of  this  thesis,  eleven  broadband  SUS  shot  transmissions 
recorded  at  the  NW  VLA  were  selected.  All  of  these  were  dropped  along  a  single  track 
from  a  point  at  latitude  40°16.2’N,  longitude  71°10.2’W  to  a  point  at  latitude  40°16.8’N, 
longitude  70°57.0’W.  The  sampling  rate  was  3906.25  Hz,  the  duration  of  one  record  was 
8.3886  sec  and  the  number  of  samples/record  was  32768.  The  arrangement  and  the  location 
of  flie  acoustical  elements  used  for  this  study  are  shown  in  Figure  1 .  This  region  of  the 
experimental  site  has  a  nominal  water  depth  of  ~  90  m  with  little  variation  (<10m)  from  the 
VLA  to  the  positions  of  the  selected  SUS  drops.  The  receiving  VLA  had  16  hydrophones 
of  which  only  the  upper  8  successfiilly  recorded  the  SUS  transmissions  with  an  upper  depth 
of  30.5  m  and  inter-element  spacing  of  3.5m.  A  typical  soimd  speed  profile  measvued  near 
the  trial  site  and  used  in  this  analysis  is  shown  in  Figure  2. 
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Acoustical 

Latitude 

Longitude 

Distance  to 

Element 

VLA 

40°22.4’ 

71°13.5’ 

VLA  (km) 

SUS-S, 

40'’16.2' 

7ri0.2' 

13.00 

SUS-S2 

40°16.2' 

71°08.9’ 

14.30 

SUS-S3 

40°16.2’ 

71°07.6’ 

15.85 

SUS-S4 

40°16.2' 

71°06.3' 

17.60 

SUS-S5 

40°16.2' 

71°05.0' 

19.50 

sus-s« 

40°16.2’ 

7r03.7’ 

21.50 

SUS-S7 

40°16.2' 

71°02.4’ 

23.50 

SUS-Sg 

40'’16.2' 

71.01.1' 

25.70 

SUS-S9 

40°16.2' 

70°59.8' 

27.80 

SUS-S, 0 

40°16.2' 

70°58.5' 

30.00 

SUS-S,, 

40°16.2' 

70°57.0' 

32.60 

Figure  1 .  Arrangement  and  Location  of  Acoustic  Elements  Used  at  the  Experimental  Site. 


Depth  (m) 
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in.  ANALYSIS  METHODS 


A.  INVERSION  OF  BOTTOM  SOUND  SPEED 

In  determining  the  bottom  sound  speed  a  method  was  performed  to  determine  the 
critical  angle  of  the  seabed  from  the  vertical  directionality  of  the  ambient  noise  in  the  water 
colimm  (Buckingham  and  Jones,  1987).  The  properties  of  the  ambient  noise  that  we  are 
specifically  concerned  with  are  the  vertical  directionality  and  the  spatial  coherence  of  the 
fluctuations  of  two  sensors  vertically  displaced  in  the  noise  field.  These  two  quantities  are 
related  through  a  Fourier  transform  relationship,  which  is  given  by 

r(Q)  =  exp[zQ  cos(^)]  sm(0)d0  (1) 


where  0  is  the  polar  angle  measured  down  from  the  zenith  and  F(0)  is  a  dimensionless 
directional  density  fimction  which  represents  the  distribution  of  the  noise  power  in  the 
vertical.  F(0)  can  be  normalized  according  to  the  condition 

^j[F(0)sin(0)d0  =  l.  (2) 

This  indicates  that  in  the  limit  as  Q  ->  0  the  value  of  the  coherence  function  goes  to  unity. 

Buckingham  and  Jones  provided  a  simple  example  for  isotropic  noise  in  which  F(0) 
is  independent  of  0.  According  to  Eqs.  (1)  and  (2)  the  coherence  function  in  this  case  is 
shown  to  be 


sin(Q) 

Q 


where  Q  is  a  dimensionless  angular  frequency  defined  as 


c 


W 


5 


(3) 

(4) 


(0  is  the  angular  frequency,  is  the  sormd  speed  in  the  water  column,  and  /  is  the 
separation  between  the  two  sensor  positions.  This  example  is  shown  in  Figure  3. 


7 


Figure  3.  Example  of  Ideal  Vertical  Coherence  Function  Versus  Frequency. 

This  function  predicts  that  the  coherence  of  the  noise  at  two  sensors  is  real  due  to 
symmetry  of  the  noise  field  about  the  horizontal.  For  an  asymmetrical  noise  distribution  in 
the  vertical,  the  coherence  function  is  complex  with  a  large  imaginary  component.  A 
S)mimetrical  noise  distribution  is  also  predicted  by  Buckingham  and  Jones’  theoretical 
model  of  shallow  water  noise,  where  the  coherence  function  for  shallow  water  noise  is  an 
oscillatory  function  that  decays  with  increasing  Q .  At  a  given  frequency,  the  directional 
density  function  of  the  noise  predicted  by  the  model  is  closely  approximated  by 


where  a,  is  independent  of  0.  Then,  the  noise  propagates  in  directions  close  to  the 

horizontal,  with  no  acoustic  energy  traveling  with  grazing  angles  greater  than  (the  critical 
grazing  angle).  The  contribution  of  local,  wind-driven  sources  to  flie  overall  noise  field  is 
approximately  isotropic,  that  is,  the  directional  density  function  of  the  continuous 
component  of  the  noise  can  be  expressed  as 

F2(0)  =  a2,  for O<0<7r,  (6) 
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where  is  independent  of  0.  By  superposing  Eqs.  (5)  and  (6)  we  have  a  representative 
directional  density  function  for  the  total  noise  field  in  shallow  water; 


for 

V 

otherwise. 


( Tl  ^ 

<0< 

-  +  0C 

(7) 


Figure  4  shows  a  sketch  of  the  vertical  distribution  expected  in  a  shallow  water 
environment.  Through  the  normalization  condition  in  Eq.  (2),  the  critical  grazing  angle  can 
be  expressed  in  terms  of  the  other  two  parameters  as 

“2 

and  the  model  coherence  function  of  the  noise  represented  by  the  distribution  in  Eq.  (7)  is 
found  from  Eq.  (1)  to  be 


r,.(e)  =  ^sin(Q)  +  ^^^^sin[Qsin(eJ].  (9) 

Using  this  model,  the  spatial  coherence  of  two  vertically  displaced  sensors  of  the 
NW  VLA  was  computed  fi'om  extracted  ambient  noise  data.  By  spatial  coherence  is  meant 


Figure  4.  Model  of  the  Vertical  Distribution  of  Ambient  Noise  in  a  Pekeris  Waveguide. 
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the  time-averaged  crosscorrelation  coefficient  of  the  noise  observed  by  hydrophones 
separated  vertically  in  die  water  column,  i.e.. 


where  i^(co)  represents  the  cross  spectral  density  as  a  function  of  frequency  for  two 
separated  phones  and  /^(o)  and /^(co)  represent  the  spectral  density  for  each  one  of  the 
separated  phones.  To  search  the  two  unknown  parameters  a^axvd  G^in  Eq.  (9),  an 
optimization  procedure  was  implemented  by  finding  a  least  squares  fit  between  the  real  part 
of  die  measured  coherence  function,  ,  and  the  theoretical  coherence  function,  By 
defining  the  ambiguity  function 

®=Slr„(aa,.e.)-Reai{r^(Q)}p ,  (ii) 

tt=} 

values  for  and  0^  are  obtained  when  O  is  a  miniTmim 
B.  INVERSION  OF  BOTTOM  DENSITY 

To  determine  the  bottom  density,  the  concept  of  an  effective  depth  for  an  ideal 
wavegmde  was  employed  (Chapman  and  Ward,  1989).  By  considering  the  phase  change  of 
a  plane  wave  reflected  at  a  fluid/fluid  boundary,  Weston  introduced  the  concept  of  the 
effective  boundary  depth,  arguing  that  one  could  view  the  reflection  as  taking  place  at  an 
imaginary  pressure  -  release  boundary  located  at  a  specific  depth  below  the  tme  boundary. 
He  proposed  that  the  normal  modes  of  the  Pekeris  model  should  be  similar  to  those  of  an 
ideal  pressure  -  release  waveguide  whose  depth  was  equal  to  the  tme  water  depth  plus  the 
effective  boundary  depth.  The  normal  mode  wavenumbers  for  the  Pekeris  model  can  then 
be  estimated  easily  without  numerical  iteration. 

The  Pekeris  model  consists  of  a  uniform  water  layer  of  constant  depth  bounded 
above  by  a  pressure  -  release  surface  and  below  by  a  uniform  fluid  half  space  having 
different  density  and  soxmd  speed.  Consider  a  plane  wave  of  frequency /in  a  fluid  medium 
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of  sound  speed  incident  upon  a  plane  boundary  with  another  acoustic  medium  at  angle  0 

relative  to  grazing,  as  depicted  in  Figure  5. 


The  ray  segments  s  and  s  ’  are  given  by 


s  = 


AH 

sinO 


and 


e' 


2A/fcos"e 
sin0 


(12) 


(13) 


and  the  phase  of  the  reflected  plane  wave  at  points  A  and  B  is  given  by 

+  ^  (14) 

and 

^,=2k^s+^',  (15) 

CO 

in  which  k^=  — .  Requiring  <|)^  =  (l)^  the  effective  depth  of  the  boundary  is  defined  by 


2A:^sin0 


(15) 
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For  an  isospeed  water  eolumn  with  sound  speed  and  density  overlying  a 
homogeneous  solid  medium  with  compressional  soxmd  speed  c^,  shear  speed  ,  and 
density  Chapman  and  Ward  (1984),  have  shown  (assuming  <  c^and  bottom 
attenuation  is  negligible)  that  the  effective  depth  is 

H,  =  H+ 

where  iy  is  the  water  depth  and  0^  is  the  grazing  critical  angle.  If  then  Eq.  (17)  is 

reduced  to 


1-2 


-\2 


_g_ 

k... 


sin0. , 


(17) 


H,^H+ 


g 

^wsm0,  ■ 


(18) 


The  boundary  conditions  for  this  ideal  waveguide  require  the  pressure  to  vanish  at 
the  upper  and  lower  interfaces.  Correspondingly,  eigenfunctions,  \|/(z),  must  satisfy 
(Medwin  and  Clay,  1998) 


v(2)Uo=  ^  and  xk(^)L^  =  0 . 

Solutions  of  the  eigenfunctions  for  these  boundary  conditions  are 


v|/(z)=sin 


^  nnz^ 


(19) 


(20) 


and  the  vertical  component  of  the  wavenumber,  y,  must  then  satisfy  the  quantization 
condition 


Yn=^,  (21) 

where  n  is  the  integer  and  designates  the  mode  number.  The  horizontal  values  of  the 
wavenumber  are  then  given  by 

^«=fe-y'r.  (22) 
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and  the  requirement  that  Ar„be  real  gives  the  mode  cut-off.  By  substituting  y„  from  Eq. 
(21)  and  k^=  —  ,  K„  can  be  written  as 


•  n2 

..  ® 


nn 


(23) 


Furthermore,  the  group  velocity  of  the  n*  mode  is  given  by 

da 


= 


dK' 


(24) 


therefore, 


M. 


=£i 

(0 


00 


"w  J 


V 


nn 


kH.j 


=  c.Jl- 


nnc^ 

(oH„ 


(25) 


The  group  velocity  u„  and  the  mode  travel  times  are  related  by 

R 


(26) 


where  R  is  the  horizontal  distance  between  source  and  receiver  and  t„.  is  the  travel  time  of 
the  n*  mode.  Combining  Eqs.  (25)  in  (26)  the  mode  travel  times  can  be  written 

1 


t  =— =— 
M. 


"  nc 


(27) 


1- 


From  the  definition  of  effective  depth,  this  becomes 


c„. 


1- 


-|2 


n7tsin0„ 


^  InfH  sin0,  ^ 


-1/2 


(28) 


which  is  the  final  expression  for  the  mode  travel  times  .  Since  we  know  the  real  depth  H, 
and  and  0^  were  determined  from  the  previous  analysis,  then  we  can  invert  for  g  by 
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measuring  the  mode  travel  times  over  a  range  of  frequencies  for  different  SUS 
transmissions.  Assuming  « 1,  then  g  «  p^. 

C.  INVERSION  OF  BOTTOM  ATTENUATION 

Mode  attenuation  as  a  function  of  range  can  be  estimated  from  several  computed 
SUS  spectrograms.  The  mode  attenuation  per  km,  a„ ,  satisfies  the  relation 

ATL„=a„AR,  (29) 

where  and  AR  are  the  mode  transmission  loss  due  to  attenuation  and  range  between 
SUS  drops,  respectively.  Neglecting  volume  attenuation  in  the  water  column,  the  decrease 
in  mode  amplitude  is  expected  to  be  caused  only  by  bottom  attenuation  and  cylindrical 
spreading  =  lOlog/?).  Since  mode  n  can  be  thought  of  as  due  to  rays  propagating  at 

3ngle  ,  the  mode  attenuation  can  be  considered  the  result  of  multiple  bottom  reflection 
losses. 


The  number  of  bottom  boxmces  over  1  km  of  horizontal  propagation  can  be  found 
by  considering  the  angle  of  elevation  or  depression  0  of  the  local  direction  of  propagation 
of  the  waveform.  Referring  to  Figure  6,  the  angle  0  and  the  diagonal  distance  for  one 
bounce  Z,  are  related  by 


sin0  = 


2H 


(30) 


where  H  is  the  water  depth.  Defining  the  total  diagonal  distance  L  corresponding  to  i?  =  1 
km  of  horizontal  propagation  by 


L  = 


R 


COS0 


(31) 


then  L  and  Z,  can  be  related  in  order  to  define  the  number  of  bottom  bounces  per  km  as 


Z  _  /?tan0 
X~  IH 


(32) 


The  attenuation  for  mode  n  can  then  be  expected  to  be  due  to  n{b„  )  bottom  boxmces  pw 
km  producing 

a„=v(0„)a,,(0j,  (33) 
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where  0„  is  the  corresponding  propagation  angle  for  mode  n  and  ctjj(0„)is  the  bottom  loss 
per  bounce  for  a  single  reflection  at  this  angle. 


The  reflection  coefficient  at  a  fluid-flmd  interface  for  a  plane  wave  traveling  at 
angle  0  is  defined  by 

^  p,c,smQ+ip^b2  ^ 

PiCjSin0-/p>2  ’ 

where  and  are  the  soimd  speed  and  density,  respectively,  in  the  upper  medium,  and 
Pj^  are  the  soimd  speed  and  density,  respectively,  in  the  lower  medium,  and 


^2  = 


r  -  ^ 


-11/2 


cos^  0-1 


(35) 
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If  there  is  no  attenuation,  then  |R|  =  1  for  incidence  beyond  the  critical  angle.  However,  if 
the  bottom  is  lossy  then 

M  =  (36) 

and  the  bottom  loss  per  bounce  can  be  defined  as 

aw(0j=-2Olog(|i?(0„)|)s2Os(0„).  (37) 

Combining  the  equations  above,  the  reduction  in  the  reflection  coefficient  amplitude  can  be 
determined  by 

"v  1  ct  1  A7Z,„ 

"  '’2OiV(0j~^A/?Ar(0j' 

Bottom  attenuation  can  be  introduced  by  defining  an  imaginary  component  of  the 
bottom  sound  speed. 


(39) 

such  that  the  wavemunber  in  the  lower  medium  becomes 

K~  ^  7  (40) 

where 


(o6 


Since  the  solution  in  the  lower  medium  goes  like 


(41) 


then  Oe  is  found  to  produce  the  exponential  decay  in  range, 
loss  according  to 


(42) 

This  then  causes  transmission 


^’^^atten  =  8.686a^M = aM ,  (43) 

where  a  is  the  bottom  attenuation  in  units  dB/m,  and  can  be  written  as 

a  =  8.686—.  (44) 

In  order  to  estimate  the  bottom  attenuation  (dB/m),  a  simple  matching  algorithm 
may  be  employed  which  searches  over  values  of  5  (the  complex  bottom  sound  speed)  until 
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the  reduction  in  the  reflection  coefficient  amplitude  agrees  with  the  observed  value  of  £(0  J 
defined  in  Eq.  37.  Once  values  for  8  are  obtained  (for  various  frequencies  and  mode  angles 
0  then  Eq.  43  may  be  employed  to  produce  the  final  estimate  for  the  bottom  attenuation 
a. 
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IV.  DATA  ANALYSIS  AND  RESULTS 


A.  INVERSION  OF  BOTTOM  SOUND  SPEED 

As  described  in  Section  HI.  A.,  the  inversion  for  bottom  sound  speed  is  based  on  the 
evaluation  of  the  vertical  coherence  of  the  ambient  noise  field.  For  this  analysis,  the 
ambient  noise  was  extracted  from  segments  of  the  SUS  data  sets  between  actual  SUS 
signals  (and  allowing  for  reverberation  to  fade).  The  spectrum  of  the  ambient  noise 
computed  showed  little  variation  from  low  frequency  up  to  600  Hz  indicting  no  strong 
unexpected  noise  sources  were  present.  Also,  the  imaginary  component  of  the  coherence 
between  the  two  hydrophones  was  found  to  be  small  as  required.  This  indicates  that  the 
ambient  noise  field  at  the  trial  site  for  frequencies  below  600  Hz  has  a  symmetrical 
distribution  in  the  vertical  direction  and  the  data  is  suitable  for  inversion  by  the  method 
described  previously. 

It  should  be  noted  that  the  inversion  method  described  in  Section  HI.  A.  was 
originally  defined  for  an  isospeed  water  column  overlying  a  homogenous  bottom.  As  is 
obvious  from  the  local  sound  speed  data  (refer  to  Figure  2),  there  is  significant  variation  in 
sound  speed  over  deptii.  However,  if  the  sound  speed  change  is  not  too  large  over  the 
vertical  separation  of  the  phones,  the  inversion  method  should  still  be  applicable  by  simply 
adjusting  the  local  angles  to  the  ^propriate  values  at  tiie  bottom  interface  using  Snell’s  law. 

Since  the  data  from  the  upper  two  phones  of  the  VLA  was  analyzed,  a  local  soimd 
speed  of  c^  «  1485  m/s  was  used  in  the  calculation  of  the  ambiguity  fimction  (refer  to  Eq. 
11). 

The  fit  between  the  theoretical  coherence  obtained  from  the  model  and  the  real 
component  of  the  observed  coherence  is  shown  in  Figure  7.  By  computing  the  ambiguity 
surface  for  the  optimization  of  two  parameters  (Oj  and  0^ )  an  optimal  value  was  formd  at 

0^«18°  in  the  water  column.  Figure  8  shows  the  ambiguity  surface  for  the  optimization 
with  two  parameters,  aj  and  0,.  whose  minimum  corresponds  to  aj  *  018  and  mean  critical 
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angle  0^«18°±r  for  4  different  SUS  explosive  data.  A  grazing  eritical  angle  of 
approximately  19°±r  near  tiie  bottom  was  determined  from  Snell’s  law.  The  soimd  speed 
close  to  the  bottom  was  1476  m/s,  and  therefore  the  bottom  sound  speed  Cj^l 563+10  m/s. 
This  value  was  compared  to  values  obtained  by  a  sophisticated  gaietic  algorithm  technique 
employed  at  the  University  of  Rhode  Island  (Miller  and  Potty,  1998).  They  found  values  in 
the  upper  few  meters  of  the  sediment  in  agreement  with  the  result  reported  here. 

B.  INVERSION  OF  BOTTOM  DENSITY 

To  invert  for  bottom  density,  the  method  of  ideal  waveguide  effective  depth  was 
introduced.  However,  this  method  assumed  a  constant  sound  speed  in  the  water  column.  In 
fact,  this  was  fundamental  in  the  definition  of  the  ideal  waveguide  eigenfunctions.  For  that 
reason,  this  method  may  not  be  expected  to  produce  very  good  results.  We  will  attempt  to 
invert  for  density  anyway  by  assuming  the  water  column  is  isospeed  with  nominal  sound 
speed  »  1485  m/s.  Consistent  with  this  value,  we  shall  assume  the  critical  angle  0, »  18°. 
In  order  to  obtain  separated  mode  travel  times,  spectrograms  were  computed  for  various 
SUS  explosive  data.  The  theoretical  mode  travel  times,  for  the  ideal  waveguide  were  tiien 
computed  as  a  function  of  range  R  between  source  (SUS)  and  receiver  (VLA),  water  depth 
H,  water  sound  speed  and  critical  angle  0^.  A  comparison  between  theoretical  and 
measured  travel  times  for  all  11  SUS  data  sets  produced  an  average  value  for  the  density 
ratio  of  1.6.  However,  the  changes  in  the  theoretical  curves  for  various  density  ratios  was 
not  large  whereas  the'  mode  arrival  structures  were  fairly  ambiguous.  Thus  the  uncertainty 
is  very  large.  It  is  also  possible  that  the  tiieoretical  mode  travel  times  have  significant  errors 
due  to  the  mismatch  in  sound  speed  profile.  However,  the  estimated  value  of  «  1.6  g/cm^ 
from  the  inversion  method  seems  quite  realistic.  The  results  supplied  by  Miller  and  Potty 
suggest  a  more  reasonable  estimate  for  the  bottom  density  is  about  1.8  g/cml  Figures  9a 
and  9b  show  the  mode  travel  times  for  a  phone  from  a  SUS  explosive  data  and  those  for  an 
ideal  waveguide. 
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C.  INVERSION  OF  BOTTOM  ATTENUATION 

A  method  was  developed  in  Section  in.C.  to  invert  for  bottom  attenuation  using 
measures  of  modal  transmission  loss  as  a  function  of  range.  In  ordCT  for  this  method  to 
work  successfully,  the  decrease  in  amplitude  due  to  bottom  attenuation  must  be  observable 
over  the  transmission  loss  due  to  cylindrical  spreading.  For  example,  as  the  range  between 
the  SUS  drop  location  and  the  VLA  increases  from  15  km  to  30  km,  the  change  in  TL  due 
to  cylindrical  spreading  is  approximately  3  dB.  This  equates  to  about  0.2  dB/km.  In  order 
for  bottom  attenuation  to  be  distinguishable,  Eq.  (38)  suggests  that  the  quantity  20N  8 
should  be  of  this  same  order  of  magnitude  or  larger. 

A  simple  analysis  assuming  bottom  attenuation  increase  linearly  with  frequency 
shows  that  8  increases  with  increasing  frequency.  Specific  values  can  be  obtained  using 
Eqs.  (34)  -  (44)  with  realistic  values  of  bottom  attenuation  a.  However,  as  frequency 
increases,  the  angle  of  propagation  for  a  specific  mode  decreases  thereby  decreasing  N,  the 
number  of  bottom  interactions  per  km. 

If  one  assumes  a  reasonable  value  of  a  ~  0.05  dB/km/Hz,  a  100  Hz  signal  suffers  a 
0.5  dB/km  loss  due  to  attenuation  for  angles  of  propagation  of  about  12°.  At  this  frequency 
in  this  waveguide,  that  corresponds  roughly  to  Mode  2.  Unfortunately,  this  mode  is  not 
clearly  distinguishable  at  this  frequency  from  the  available  SUS  data.  The  lowest  modes 
can  be  isolated  at  lower  frequencies,  but  then  the  attenuation  is  negligible.  At  higher 
frequencies,  the  propagation  angles  for  these  low  modes  is  near  grazing,  and  again  the 
influence  of  attenuation  is  small. 

The  result  of  this  analysis,  unfortunately,  shows  that  this  data  does  not  avail  itself  to 
inversion  for  bottom  attenuation  using  the  method  developed  previously.  Other  processing 
methods  may  be  able  to  extract  the  necessary  information  in  order  to  utilize  this  method,  but 
will  not  be  explored  further  here. 
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Figure  7.  Coherence  of  the  Measured  Ambient  Noise  between  Two  Hydrophones;  Real 
Component  (Dashed  Curve)  and  hnaginaiy  Component  (Dotted  Curve).  A  Good  Fit 
Between  the  Theoretical  Coherence  Obtained  from  the  Model  (Solid  Line)  and  the  Real  of 

the  Measured  Coherence  is  Observed. 


Figures.  Ambiguity  Surface  for  the  Optimization  with  Two  Parameters,  and  G^.  The 
Minimum  is  Marked  by  X  which  Corresponds  to  «  0. 1 8  and  0^ »  1 8°  for  4  Different  SUS 

Explosive  Data. 
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Figure  9a.  Mode  Travel  Times  for  a  Phone  from  a  SUS  Explosive  Data, 
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Figure  9b.  Mode  Travel  Times  for  a  Phone  from  an  Ideal  Waveguide. 


24 


V.  CONCLUSION 


This  thesis  incorporated  several  methods  which  may  be  applied  directly  to  measured 
data  to  provide  a  relatively  simple  means  of  obtaining  reasonable  estimates  of  ocean  bottom 
parameters  from  minimal  information.  The  analysis  of  the  data  from  the  Mid-Atlantic 
Bight  experiment  recorded  on  the  NW  VLA  produced  the  following  results: 

1.  Ambient  noise  was  used  to  measure  both  the  vertical  directionality  and  spatial 
coherence  of  the  noise  field.  By  examining  the  noise  data,  we  have  been  able  to 
estimate  0^  and  hence  the  sound  speed  c^,  in  the  sediment.  The  bottom  sound 
speed  deduced  was  approximately  1563±10  m/s. 

2.  The  method  of  effective  depth  was  employed  to  estimate  the  bottom  density  p^. 
Although  this  method  was  derived  for  an  isospeed  water  column  environment,  it 
is  expected  to  produce  reasonable  estimates  for  most  shallow  water 
environments.  By  comparing  mode  travel  times  extracted  from  SUS  data 
spectrograms  and  comparing  these  with  values  for  an  ideal  waveguide  with  an 
effective  depth,  the  bottom  density  was  estimated  to  be  «  1.6  g/cm^.  The 
uncertainty  in  this  value  is  very  large,  however,  partly  due  to  the  difficulty  in 
distinguishing  mode  arrivals  clearly  in  the  spectrograms. 

3.  In  order  to  distinguish  bottom  attenuation  from  simple  cylindrical  spreading,  it 
was  determined  that  values  for  mode  transmission  loss  were  required  for  high 
modes  at  high  frequencies.  Unfortunately,  the  SUS  spectrogram  data  could  not 
adequately  provide  this  information,  and  no  estimate  of  bottom  attenuation  was 
produced.  It  may  be  possible  to  use  the  method  outlined  if  a  modal 
decomposition  of  the  data  on  the  VLA  could  be  performed.  However,  due  to  the 
relatively  short  aperture  of  the  array  (compared  to  (he  waveguide  depth),  this 
seems  unlikely.  Furthermore,  the  source  level  of  different  SUS  detonations  may 
not  be  repeatable  enough  to  use  this  data  for  accurate  estimations  of  bottom 
attenuation. 
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Comparing  the  inversion  results  for  bottom  sound  speed  and  density  with  values 
obtained  from  genetic  algorithm  techniques  provided  by  the  University  of  Rhode  Island, 
good  agreement  was  found  particularly  for  the  sound  speed.  However,  the  analysis  of  the 
SUS  data  proved  more  difficult  than  originally  anticipated.  Much  of  this  has  to  do  with  the 
imcertainty  in  absolute  drop  location  and  repeatable  source  signals.  Another  part  of  the 
problem  was  the  reliance  on  low  mode  information  which  is  difficult  to  characterize  for 
SUS  data  measured  on  only  eight  hydrophones  with  total  aperture  less  than  half  the  width 
of  the  waveguide.  It  is  likely  that  this  analysis  would  be  more  successful  with  different 
more  well-defined  sources  and  a  full  water  column  array. 
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